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Abstract. For all hyperbolic polynomials we proved in [11] a Lipschitz esti- 
mate of Jacobi matrices built by orthogonalizing polynomials with respect to 
measures in the orbit of classical Perron-Frobenius-Ruelle operators associated 
to hyperbolic polynomial dynamics (with real Julia set). Here we prove that 
for all sufficiently hyperbolic polynomials this estimate becomes exponentially 
better when the dimension of the Jacobi matrix grows. In fact, our main result 
asserts that a certain natural non-linear operator on Jacobi matrices built by 
a hyperbolic polynomial with real Julia set is a contraction in operator norm 
if the polynomial is sufficiently hyperbolic. This allows us to get for such poly- 
nomials the solution of a problem of Bellissard, in other words, to prove the 
limit periodicity of the limit Jacobi matrix. This fact does not require the 
iteration of the same fixed polynomial, and therefore it gives a wide class of 
limit periodic Jacobi matrices with singular continuous spectrum. 



1. Introduction 

Let T be an expanding polynomial with the real Julia set Julia(T), degT = d. 
We recall that Julia(T) is a nonempty compact set of points which do not go to 
infinity under forward iterations of T . Under the normalization 

T" 1 : [-1,1] -> [-1,1]; ±leJulia(T) (1.1) 

such a polynomial is well-defined by the position of its critical values 

CV(T) := {U = T(a) : T '(a) = 0, a > Cj for i > j}. 

Expanding, or hyperbolic polynomials are those, for which 

Ci i Julia(T),Vi, 

which is the same as to say that CV(T) n Julia(T) — (just use the fact that 
Julia(T) is invariant under taking full preimage T _1 ). The term "expanding" is 
deserved because for expanding polynomials one has the following inequality 

3Q > 1, |(T n )'(ir)| > cQ n ,Vx e Julia(T) . (1.2) 

Here and in everything what follows T n means n-th iteration of T, T n — T0T0....T. 
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Let us mention that for T with a real Julia set one has |T(ci)| > 1 since all 
solutions of T(x) — ±1 should be real. 

We will need to consider the notion of "sufficiently expanding" ("sufficiently 
hyperbolic") polynomials. As we saw, the expanding property is the same (in our 
normalization) as dist(CV(T), [—1, 1]) > 0. The polynomial T with normalization 
(1.1) will be called sufficiently hyperbolic (or sufficiently expanding) if 

dist (CV(T), [-1,1]) > A, (1.3) 

where A is a large absolute constant to be specified later (but A = 10 will work). 
Notice that the definition of sufficient hyperbolicity does not involve the degree 
of T. In particular, T and any of its iterative powers T 2 ,T 3 ,... are sufficiently 
hyperbolic simultaneously. 

A Jacobi matrix J : l 2 (Z) — ► Z 2 (Z) is called almost periodic if the family 

{S k JS k }kez, 

where S is the shift operator in ? 2 (Z), S\m) = \m + 1), is a precompactum in the 
operator topology. 

Example. Let G be a compact abelian group, p(a), q{a) be continuous functions on 
G, p(a) > 0. Then J(a) with the coefficient sequences {p(a + fc/i)}fc, {q{a + fc/i)}fc, 
/j, € G, is almost periodic. 

Let us show that in fact this is a general form of almost periodic Jacobi matrices. 
For a given almost periodic J define the metric on Z by 

pj(k) :=\\S- k JS k - J\\. 

Evidently pj(fc + to) < pj(k) + pj(m). Then J = J(0), where G — Ij, Ij is the 
closure of Z with respect to pj, and p = 1 e Ij. 

Recall that for a given system of integers {dk\kLi one can define the set 

I = lim{Z/rfi...4^}, (1.4) 

that is a € I means that a is a sequence {ao, ai, Q2, •■•} such that 

at G Z/di...dfe + iZ and afc|modrfi...rffe = «fe-i- 

The addition in I is defined as addition in the i-th entry. The metric dist(a, (3) = n l , 
where k 6 (0,1), I is the smallest integer such that a\ ^ Pi, makes I a compact 
abelian group. In particular, if p is a prime number and dk = p we get the ring of 
p-adic integers, I = Z p . 

In this work we build a certain machinery that allows to construct almost periodic 
Jacobi matrices with singularly continuous spectrum such that Ij = I. 

The key element of the construction is the following 

Theorem 1.1. Let J be a Jacobi matrix with the spectrum on [—1,1]. Then the 
following Renormalization Equation has a solution J = J(J) = J(J;T) with the 
spectrum on T _1 ([— 1, 1]): 

V*(z - jy'V = (T(z) - J)- 1 T'(z)/d, (1.5) 

where V\k) = \dk). Moreover, i/min^ |tj| > 10 then 

\\J(Jl) - J(J2)\\ < K\\J, - J 2 \\. 

with an absolute constant k < 1 (does not depend on T also). 

This theorem, for example, will result in the following consequence: 



LIMIT PERIODIC JACOBI MATRICES 



3 



Theorem 1.2. Let T be sufficiently hyperbolic in the sense of (1.3). Let be the 
Jacobi matrix obtained by orthogonalizing polynomials with respect to the balanced 
(equilibrium) measure uj on the Julia set of T . Then J w is a limit periodic matrix. 
In other words, the sequences that give the diagonal and the below (above) diagonal 
entries are uniform limits of periodic sequences. 

Remarks. 1) A known problem (due to Bcllissard) asks to prove this statement for 
all hyperbolic polynomials normalized as in (1.1). Here we do it for all sufficiently 
hyperbolic polynomials (recall that we always tacitly assume that Julia(T) is on the 
real line). Our result about sufficiently hyperbolic polynomials explains, in a sense, 
the earlier results in [2], [3], where it has been shown that polynomials T(z) := 
e~ n T n (ez) generate limit periodic Jacobi matrices if e is sufficiently small, here T n 
is the n-th Tchebyshcf polynomial. Smallncss of e obviously makes T = e~ n T n (sz) 
sufficiently hyperbolic in our sense (1.3). 

2) In the thesis of Herndon [7] Theorem 1.2 is proved by another method. We 
regret that it has not been published, that might have clarified the proof, which 
seems to be quite involved. 

3) One can wonder after analyzing the results of [2], [3] and the present paper, 
that may be there is a threshold of hyperbolicity: before it J w is not limit periodic, 
and after it it is limit periodic. However, we do not believe in this sort of behavior, 
but at this stage we cannot prove the conjecture of Bcllissard for all hyperbolic 
polynomial with real Julia set. 

4) Let us mention that, in fact, (1.5) has 2 d ~ 1 solutions such that the spectrum 
of J is on T _1 ([— 1, 1]). Here we use only one of them. 

We note that the real output of Theorem 1.1 is much wider than Theorem 1.2. 
ft shows that 

a) roughly speaking, constructing in a regular iterative way a Cantor set E, 
E C • • • C E n+ \ C E n . . . , that may support the spectrum of a limit- 
periodic Jacobi matrix it is enough to follow the strategy: on each step the 
approximating set E n should have a form of an inverse polynomial image, 
i.e.: 

E„ = f/"^— 1, 1], U„ is a polynomial, 

or, what is the same, E n should be the spectrum of a periodic Jacobi matrix; 

b) the above statement becomes a theorem if on each step we remove from the 
previous set a sufficiently large part (using sufficiently expanding polyno- 
mials), i.e.: if T\,T2-.., is a sequence of polynomials with sufficiently large 
critical values, then E n := U~ x [— 1, 1], with U n = T n o ... oT 2 o T\; 

c) the set E, that was constructed in the above described way, is the spectral 
set of infinitely many (uncountable set) of different limit periodic Jacobi 
matrices, that means that each of the matrices does not belong to the hull 
of another one (can't be obtained as a limit of shifts). The problem: to 
describe the set of all limit periodic Jacobi matrices with spectrum E or 
certain subclasses (or at least to try to understand how these sets look like), 
is a very interesting and challenging problem. 
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Let us outline a proof of claim b). First, we point out the following two properties 
of the function J(J; T) in Theorem 1.1. Due to the commutant relation VS = S V 
one gets J{S- m JS m ) = S- dm J(J)S dm . The second property is 

J(J(J;T 2 );Ti) - J(J;T 2 oTi), 

that is, the chain rule holds. 

Now, we produce a limit periodic Jacobi matrix with singularly continuous 
spectrum and such that Ij = I. For the chosen system of polynomials Ti,T2..., 
deglfc = rffe, with sufficiently large critical values, define J m — J( J; T m o ...oT 2 oTi). 
By Theorem 1.1, the limit J = lim m ^oo J m exists and does not depend on J. More- 
over, 

Vj, \\J - S- dl - dl3 JS dl - d ' 3 \\ < \\J- S-3JS 3 \\k 1 <2k 1 . 

That is pj defines on Z the standard p-adic topology in this case. This proves that 
J is a limit periodic matrix, in particular, it is almost periodic. 

Notice that for the case T\ = T 2 = ... = T m =: T, we get the limit periodic 
matrix with the spectrum on Julia(T). 



2. Renormalization equation 

In this section it is convenient to assume that 

T{z) = z d -qdz d - x + ... 

is a monic expanding polynomial. Under this normalization T^ 1 : [— — > [— 
for a certain £ > 0. 

Let J : l 2 (Z) —> l 2 {1) be a Jacobi matrix with the spectrum on [—£,£]. We 
describe the set of solutions of the Renormalization Equation 

V*(z-J)- 1 V = (T(z)-J)- 1 T'(z)/d, V\k) = \kd), (2.1) 

here J is a Jacobi matrix with the spectrum on T _1 ([— £, £]) that should satisfy 
(2.1). 

In what follows by l±(s) we denote the spaces which are formed by {|s + /c)} with 
k < and k > respectively, that is Z 2 (Z) = l-(s) © l\(s + 1). Correspondingly to 
these decompositions we set J±(s) — P x i ^ J\l±(s). 

Recall that a (finite or infinite) one-sided Jacobi matrix is uniquely determined 
by its so called resolvent function 

r±(z,s) = (s\(J±(s)-z)- 1 \s), (2.2) 

for which the following decomposition in the continued fraction holds true 



Q s - 



1 -L_ 



r+(z,s) = = - = ; - 2 z - : , tt - (2.3) 



z- q s +i - ••• 

Lemma 2.1. Assume that a matrix J satisfies (2.1). Then 

Psd+i-Psd+d=Ps+i and q s d = q- (2.4) 
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Further, let J 1 -^ be the s-th d x d block of J, that is 



q s d Psd+l 
Psd+1 q s d+l Vsd+2 



Psd+d-2 Qsd+d-2 Psd+d-1 
Psd+d-1 q s d+d-l 



(2-5) 



Then its resolvent function is of the form 



T'(z)/d 
T(')(z) ' 



(2.6) 



where T^ s \z) is a monic polynomial of degree d. Moreover, at the critical points 
{c : T'(c) = 0} the following recurrence relation holds 



TW(c) 



Pi 



T( s -!)(c) 



T{c) - q s . 



(2.7) 



Proof. We write the Jacobi matrix J as &dxd block matrix (each block is of infinite 
size), that is, we are just reordering the standard basis: 



Qo Pi 
Pi Qi 



P d S* 



SP d 



P2 



Pd-2 Qd-2 Pd-1 

Pd-i Qd-i 



(2.8) 



Here Pk (respectively Qk) is a diagonal matrix Pk = diag{pfe+sd}s>o and S is 
the shift operator. With respect to this reordering V* is the projection on the 
first-place block-component. 

Using this representation and the well known identity for block matrices 



'A B 


-l 


'(A - BD~ 1 C)~ 1 * 


C D 




* * 



we get 



T{z) - J 
T'(z)/d 



= z-Qo- [Pi, 



SPd] {z-JiY 



Pi 



PdS* 



(2.9) 



where Ji is the matrix that we obtain from J by deleting the first block-row and 
the first block-column in (2.8). Thus the second relation in (4.1) is already proved, 
we just compare the leading terms in the decomposition over powers of 1/z in the 
right and left hand sides and note that the third term on the right is of order 1/z. 

But the most important remark is that in (z — Ji) each block is a diagonal matrix 
(means all diagonals are the main diagonals in each block, on the contrary to J 
that contains SPd and PdS*). That's why we can easily get an inverse matrix in 
terms of the scalar orthogonal polynomials. 

Let us introduce the following notations: everything related to has super- 
script s. For instance: p^ = p s d+k, 1 < k < d, respectively and mean 
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orthonormal polynomials of the first and second kind, in particular, 

(«)/ 



(0|(z-J^)- 1 |0) = %^. 



P ( d S \z) 



(2.10) 



Let j[ s ^ denote the matrix that we obtain from (see (2.5)) by deleting the 
first row and the first column. Then, for j[ s ^ , Q d ^ are orthogonal polynomials of 

(s) 

the first kind and we denote by R d corresponding to them orthogonal polynomials 

(s) (s) 

of the second kind. Note that P d and R d are related by 



Q ( d a \z) 



(2.11) 



(s) 

In these terms the four interesting for us elements of the resolvent of J{ are: 



,(«) 



www 



WWW 



Q 



T7T 



(2.12) 



Recalling again that J\ is just a block-diagonal matrix we substitute (2.12) in 
(2.9). As result in the RHS (as well as in the LHS) we get a three-diagonal matrix. 
On the main diagonal we have 

z-Qo- Pi{(z - Ji)" 1 }!,!^ - SV d {(z - Ji^d-i^-iVdS*, 

and each entry on the diagonal, due to (2.12) and then (2.11), is 



(p[ S+1) ) 2 



R 



(s+l) 



s) 



Q d s \z) 



p. 



(M-l) 



,(«+!) 



Q d s \z) 



Comparing this with the LHS (2.9) we get 
T(z)-q s+1 P^ +1) (z) 



iPd) Q d s \ Z ) 



T'(z)/d Q d s+1 \z) 



(») 



(2.13) 



Similarly, below the main diagonal on the right we have 

SV d {(z- Jl)- 1 }^!,!^. 

So, using (2.12), we get from (2.9) 

1 



Ps+l 



T'(z)/d q^ {z) z« 



W (») 

gl -Pd 



(2.14) 



Thus the first relation in (4.1) is also proved, moreover all Q d s \z) (independently 
on s), being normalized to the leading coefficient one, coincides with T'(z)/d. 
We define T<- S \z) = z d + ... by the same normalization 

Tl'\z):=p a+1 P<f\z). (2.15) 

Then (2.10) implies (2.6). 
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Now we use the (last) well known fact on orthogonal polynomials — the Wron- 
skian identity: 

due to which, if T'(c) = then 

-P { MUc) = -Jy-- (2-16) 

p d ( c ) 

So, combining (2.16) with (2.14), we get from (2.13) the recurrence relation 

T(c)-q s+1 =T^\c) + J^. (2.17) 

Thus the lemma is completely proved. □ 

Now we are in a position to show that the Renormalization Equation has 2 d_1 
solutions. Then we show that they are the only possible solutions. This set of 
solutions we parametrize by a collections of vectors 

5 := {5 C } C1 (2.18) 

where each component 5 C can be chosen as plus or minus one. 

Theorem 2.2. Fix a vector S of the form (2.18). For a given J with the spectrum 
on [—£, £] define the Jacobi matrix J according to the following algorithm: 
For s € Z we put 

-l^. = -r_(T(c),s), ifS c = -l, (2.19) 

and 

T^(c) = -p 2 s+1 r + (T(c), S + l), ifS c = l, (2.20) 
where the functions r±{z, s) are defined by (2.2). Then define the monic polynomial 
(z) of degree d by the interpolation formula 

T^(z) = (z-q)T'(z)/d + £ _^£L T « (c) . (2.21) 

c:T'(c)=0 ^ / W 

Define the block (see (2.5)) by its resolvent function according to (2.6). Finally 
define the entry p sd+d by (4.1). 

We claim that the matrix J = J(S,J), combined with such blocks and entries 
over all s, satisfies (2.1). 

Example. The solution related to the vector 

5_ ={-!>• •-,-!}, 

that is all T^(c) arc defined by (2.19), plays the most important role in what 
follows. Precisely for this solution J (J) := J(6-,J) we prove the contractibility 
property (our main Theorem 1.1). 

Proof. First of all let us mention that for all c, \T(c)\ > £, that is T^(c) is well 
defined by either (2.19) or (2.20), moreover this value is of the same sign as T(c). 
It means that the rational function — ^ff^ is the Stieltjes transform of a positive 

measure (supported at the zeros of T^(z)), and hence there exists a unique d x d 
Jacobi matrix defined by (2.5). Note that (4.2) implies immediately that q s d = q, 
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again we look at the leading term in the decomposition of the resolvent function 
into the continued fraction. 

Now, we see that (2.14) holds, it's just a matter of the definition of . We 

have to check (2.13). Using the form of Q^\z) we have equivalently 

T(z) - q s+1 = T^(z)+p s+lP ( d s) Q { d s ! 1 (z). (2.22) 

Subtracting (z — q)T'(z) / d from both parts we arrive at the question of the identity 
of two polynomials of degree d — 2. Thus, we need to check (2.22) only at the 
critical points. Using the Wronskian identity (2.16) we get (2.17). Of course, the 
main point is that T^ s '(c), being defined by either (2.19) or (2.20), satisfy the 
recursion (2.7) (this is the continued fraction decomposition for f±(z, s), see (2.3)). 
Thus having (2.9), we proved (2.1). 

□ 

Theorem 2.3. Theorem 2.2 describes the whole set of solutions of the Renormal- 
ization Equation. 

Proof. We need to show that (2.19), (2.20) give the complete possible choice of the 
values T( s )(c), say for s = 0. Then all other values are of the same form due to 
Lemma 2.1. 

We claim that any other choice of T^(c) contradicts to the regularity of the 
resolvent of J at c. 

Using standard formulas for orthogonal polynomials for two-sided Jacobi ma- 
trices (see Appendix, Corollary 7.2) we have from the Renormalization Equation 

rZ 1 (T(z),0)=rZ 1 (z,0)r(z)/d + pjp 1 Rf\z), 
p\r+{T{z),l) =p\r + {z,l)T'{z)/d + p\p 1 Rf{z), 

where r±(z,s) are the resolvent functions of J±(s). 

In the same time both functions r + (z, 1) and rZ 1 (z, 0) can not have a pole at c 
simultaneously. This contradicts to (see (7.5)) 

(1|(J - ^)- 1 |l> = r + (z, ljrzHz, 0)(0\(J - z)- 1 ^) : (2.24) 

(0|(J — z) _1 |0) can not have zero multiplicity more than one and (1|(J — z) _1 |l) 
should be regular at c. 

By (2.11) T(°)(c) = -piPiR { p{c), so wc get (2.19) from (2.23) if rz\z,0) is 
regular and (2.20) in the case when r + (z, 1) is regular at c. □ 

Remark. Formulas (2.19), (2.20) play the main role in solving the Renormalization 
Equation. Actually, we proved or found them in Theorem 2.3 for s = 0, or, in the 
same way, for any other fixed s. However, there is another important ingredient: 
it should be also shown that if we choose this or that definition of (c) for any 
particular s we have to be consistent with this definition for all other values of s, 
that is we can not define, for instance, T^(c) by (2.19) and TW(c) by (2.20) for 
the same c. Precisely this claim is the main output of Lemma 2.1. □ 
Let us mention that the Renormalization Equation can be rewritten equivalently 
in the form of polynomials equations. 



Lemma 2.4. Equation (2.1) is equivalent to 

V*T(J) = JV\ 



(2.25) 
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V* T ^ _^ V = T'(z)/d. (2.26) 

Proof. Starting with (2.25), (2.26) we get 

(T(z) - J)V*{z - J)- X V = V*{T(z) - T{J)}{z - jy x V = T\z)/d. 
Having (2.1) we get 

V* T{Z) ~_ T J ( - J) V =T(z)V*(z - J^V - V*^jV 

=T( Z )lMK__ V *m T V (2.27) 
T{z)-J z-J v 1 

=T'{z)/d + JV*(z - J)- l V - V*^PjV. 

Since the left hand side in (2.27) is a polynomial of z we obtain two relations 

v , T{z)-T{J) v = T 

z — J 

and 

{ jv* - v*t(j)}{(z - jy^} = 0. 

Since vectors of the form (z — J)~ 1 V f ', / e I 2 , are complete in I 2 the last relation 
implies (2.26). □ 

3. Proof of the main theorem 

We start with 

Lemma 3.1. Let be the s-the block of J = J(5,J). Define the measure 
by 

T( s \z) f da^lx) 

— - = z — a 

T'(z)/d q 

that is, er( s ) is the spectral measure of the obliterated matrix normalized by 

J da^(x) = (p[ s) ) 2 . 

Then 

J Ti*){c) 2 p 2 +1 ■ ^ 
Proof. Note that is the orthonormal polynomial with respect to a^ s \ 

J(Q< k s l 1 (c)) 2 da^(c) = l. 
Using (2.16) and the normalization (2.15) we get (3.2). □ 
The proof of the theorem is based on the following well-known and simple lemma. 



f^^, (3.1) 
J x - z 
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Lemma 3.2. Assume that two non-normalized measures a and a are mutually 
absolutely continuous. Moreover, that da = f da and (1 + e) _1 < / < (1 + e). Let 
us associate with these measures Jacobi matrices J — J (a), J = J (a). Then for 
their coefficients we have 

\p s -Ps\ < e\\J\\, s > 0. 

Proof. Assume that p s > p s . Let us use an extremal property of orthogonal poly- 
nomials, 

(1 + e)pl...f s =(1 + e) J p 2 Q ...? s P* da>J {z s + ...} 2 da 

> inf / P 2 da = p 2 ...p 2 s . 
{P= Z °+...}J 

(l + e)po-Ps_i >pI-pI-i- 
P 2 S <P 2 S <(1 + efp 2 s 
0<Ps-Ps< CPs- 



Similarly 
Therefore 
and hence 



□ 

Proof of Theorem 1.1. Given Jj and J/j let us compare the blocks jj^ and Jjf of 
the matrices Jj := J(5-, Jj) and Ju := J(5-, Ju). Actually we will apply Lemma 
3.2 to non-normalized spectral measures a^ and a\ S j , see (3.1), corresponding to 

the obliterated matrices {Jj S ^)i and (Jjf)i- 

Note that both measures are supported on the critical points {c : T'(c) = 0} 
and, therefore, they are mutually absolutely continuous. Moreover, the density of 
the second measure with respect to the first one is of the form 



/(c) = /«(c) :=5#- 
T\ s \c) 

Assuming /(c) > 1 let us estimate /(c) — 1 from above. 

„ . _ , _ l/T,'"(e)-l/r';'(e) 

^ . . . (3.3, 

_ (s\(T(c) J ff ,_( 8 ))-i(J f ,_( s ) - J H ^(s))(T(c) - JiAsjT^s) 

1/T«(c) 

Since the spectrum of Jj/ ; _(s) is on [—£,£] we get, by definition (2.19), 

1 

> 



T#(c) 
and, for the same reason, 



\\(T(c) A-)- 1 ^ < |T(c)| _ r i = I or i = II. 
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Therefore 

0<f(c)-l<\\Ji-Jn\\ ^l^ r (3.4) 
Thus, by Lemma 3.2, we obtain 

\{pi)sd+k ~ (j>n)sd+k\ < x\\Ji ~ Jn\\, l<k<d-l, (3.5) 

where 

\T(c)\/Z + l 

x := max ■ 



(\T(c)\/t-iy 

We have to estimate \(pi) s d+d — {Pu)sd+d\- Note that due to (4.1) and Lemma 

3.1 

1 {Pi)( s +i)d . 1 . . 

^ \rrt m Tt — 7' 1 = 1 oxi = II. (3.6) 



{Pi)sd+l-{Pi)sd+d-l Ps+1 l T ( c )l/C-l 

Now, 

(Pl)s+1 (Pll)s+1 



{Pl)ad+d — {Pll)sd+d 



+ 



{Pl)sd+l--{Pl)sd+d-l ipil)sd+l--(Pll)sd+d-l 

(PfWl ~ (Pll)s+1 
(Pl)sd+l---(Pl)sd+d-l 

(Pll)s+1 (, {Pl)sd+l — (Pl)sd+d-l 



(pi)sd+l---(pi)sd+d-l V (pil)sd+l — {pil)sd+d-l 

Using (3.6), (3.4) and Lemma 3.2 we obtain 

\{pi)sd+d - (pn)sd+d\ <\\Ji - J 2 ||max 1 — - 

c \T{c)\/£-l 

+ m c ax |r(c)|/g-i l|Jl _ J2|k/2 (3J) 

H|Ji-J 2 ||niax |r(c)| 1 /c _ i (l + V2). 

Thus (3.5) and (3.7) show that say for min c \T(c)\/(, > 10 the renormalization is a 
contraction. □ 

4. Proof of Theorem 1.2 
Proof. For a given sufficiently hyperbolic polynomial T we define 

Jn+l = J(S-, Jn) 

starting from an arbitrary initial Jq — J with the spectrum on [—£,£]. Due to 
the contractibility of the renormalization J n converges in the operator norm to J, 
moreover J has the spectrum on Julia(T) and it is limit periodic 

\\J-S- dni JS dni \\ < 2^k", k<1. 

We claim that J is an orthogonal sum of two one-sided Jacobi matrices 

J_(-l) 

o j+(o)_ • 

That is, we claim that p = 0. Indeed, by Lemma 3.1 we have 

p(n+l) sd < jmax ^_ l |p(n) a < np{n) s , (4.2) 



J 



(4.1) 
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where p(n) s is the s-th coefficient in the matrix J n . Therefore, all sd l -th coefficients 
of J are uniformly small 

p sd ' < 2C«', 

and, in particular, po = 0. 

Thus J + := J+(0) is a one-sided Jacobi limit-periodic matrix with the spectrum 
on Julia(T), moreover, its spectral measure <r + (supported on Julia(T)) possesses 
the renormalization property 

This means that er + is an eigen-measure for the Ruelle operator L*a + = <r+, where 
the operator L acts on a continuous function / on Julia(T) by 

{Lf){x)= l - f(y)- (4-4) 

T{y)=x 

In other words <r + is the balanced measure on Julia(T). 

□ 



Note that due to the Renormalization Equation the spectral measure of J_ is 
the eigen-measure for the Ruelle operator 

(**/)(*)= E (4 - 5) 

T(y)=x 

i.e. L* 2 <J- = per-, p > 0. In the case of quadratic polynomials this fact was proved 
in [10]. 

5. The renormalization of periodic matrices 

The renormalization (2.1) acts in the most natural way on periodic Jacobi ma- 
trices. We recall some basic facts from the spectral theory of such matrices. 
The spectrum E of any periodic matrix J is an inverse polynomial image 

E = U~ 1 [-l,l] (5.1) 

the polynomial U of degree g + 1 should have all critical points {cjj} real and for 
all critical values |t7(cj/)| > 1. For simplicity we will assume |f (cj/)| > 1. Then 
the spectrum of J consists of g intervals 

E = [b ,a o ]\{U g j=1 (a j ,b j )). 

Also it would be convenient for us to normalize U by a linear change of the variable 
such that b = — 1 and a = 1. 

Having the set E of the above form fixed, let us describe the whole set of periodic 
Jacobi matrices J(E) with the given spectrum. To this end we associate with U 
the hyper-elliptic Riemann surface 

X = {Z= (z, A) : A 2 - 2U(z)X +1 = 0}. 

The involution on it we denote by r, 

rZ:=[z,^jeX. (5.2) 

The set 

X+ = {Ze X : \X(Z)\ < 1} 



LIMIT PERIODIC JACOBI MATRICES 



13 



we call the upper sheet of X. Note X+ ~ C \ E, in fact, z{Z) £ C \ E if Z £ X+. 

The following well known theorem describes J(E) in terms of real divisors on X. 
The Jacobian variety of X 7 Jac(X), is a g dimensional complex torus, Jac(X) ~ 
C 9 /L(X), where L is a lattice (that can be chosen in the form L = Z 9 + QZ 9 with 
Imfi > 0). Consider the g dimensional real subtorus consisting of divisors of the 
form 

9 

D{E) = {D = D+- D c , D+ := ]T Z { :Z,eI, z{Z t ) £ [a u 6,]}, 

»=i 

here -De is a point of normalization that we choose of the form 

9 

D c :=^2d:Ci€ X, z(d) = (c^, |A(Cj)| > 1, 

i=l 

- the collections of the points on the lower sheet with the z-coordinates at the 
critical points. (At least topologically, it is evident D{E) ~ R s /Z s ). 

Theorem 5.1. For given E of the form (5.1) there exists an one-to-one corre- 
spondence between J{E) and D(E). 

Let now U be a polynomial of the above described form, we restore the nor- 
malization T _1 : [—1,1] — * [—1,1] for the expanding polynomial, and we define 
U = U o T. Then we have a covering it of the Riemann surface X associated to U 
by the surface X associated to U : 

ir(z,X) = (T(z),\), (5.3) 

note 7r : X + — > X + . 

According to the general theory, this covering generates different natural map- 
pings [9], in particular, 

tt» : Jac(A) -» Jac(A), (5.4) 

and 

7T* : Jac(A) -» Jac(X). (5.5) 

In this section we solve equations (2.25), (2.26) using this language, see Theorem 
5.8. Note that (2.25) already guarantied that J £ J(E) implies J £ J(E). 

To continue we need to recall some special functions on hyper-elliptic Riemann 
surfaces. 

The first object is the Complex Green's function. Note that the function A in X + 
has no zeros except for infinity, where it has a zero of multiplicity g + 1 , moreover 
|A| — 1 on dX + . We define the Complex Green's function (with respect to infinity) 
by b 9+1 = A. It is not single valued in X + but it has the only simple zero at infinity. 
Note that 

G{z)=log m\' 

where G(z) = G(z, oo) is the standard Green's function for the domain C \ E. 
Generally, 

G(z,z ) = log 

\K K z ) I 

defines the Complex Green's function b ZQ with the only zero at Zq £ X+, z(Zo) = 
z a £C\E. 
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Since A o tt = A we have the relation 

boir = b d (5.6) 

The differential ^dlogfo, being restricted on dX + , is the harmonic measure du> 
of the domain C \ E with pole at infinity. 

The space L p (dX + ), in a sense, is the L p space with respect to the harmonic 
measure, but it should be mentioned that dX + = (E — iQ)\j(E+iO), i.e., an element 
/ of L p (dX + ) may have different values f(x + iO) and f(x — iO), ie£ 

Having in mind (5.6) we get 

f{Z) \ (Z)du (5.7) 
for every / £ L 1 {dX+). 

Definition 5.2. The Hardy space H 2 (X+) consists of functions / holomorphic on 
X + (or what is the same in the domain C\E) having harmonic majorant 

\f(z)\ 2 <u(z), zeC\E, (5.8) 

where u(z) is harmonic in C \ E. The norm of / is defined by 

H/ll 2 := inf ^(oo), 

u 

where u runs over all harmonic functions satisfying (5.8) 

An equivalent way to define H 2 (X + ) is to close the set of holomorphic functions 
uniformly bounded in X + with respect to the norm 

\\f\\ 2 --=! l/l 2 ^- (5-9) 

JdX + 

As it follows directly from (5.7), the covering (5.3) generates an isometrical 
enclosure 

v+ : H 2 (X) -» H 2 (X+) (5.10) 

acting in a natural way 

(v + f)(Z) = /(tt(Z)). (5.11) 

Now we have to describe the most complicated but the most important clement 
of the construction: we have to introduce a very natural orthonormal basis in 
H 2 (X + ). The multiplication operator by z, with respect to this basis, will lead 
us to Jacobi matrices, the substitution (5.11) to the isometry V and so on... This 
basis is a counterpart of the standard basis of {C"}n>o in the standard Hardy space 

# 2 (o),e = {|c|<i}. 

Note that 1 <E H 2 (X + ), moreover, 

</.!>=/(«>) 

for every / £ H 2 (X + ). Therefore the orthogonal complement to 1 consists of 
functions with /(oo) = 0. Let us give an alternative description of 

H 2 (X + ) = {/ e H 2 (X + ) : /(oo) - 0}. 

Any function from H 2 (X + ), having zero at infinity, is the form / = bf. However 
b is not single- valued, thus so is /. We need to generalize slightly Definition 5.2. 
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Definition 5.3. Let F = T(E) be the fundamental group of the domain C\E. Let 
a be an element of the dual group of characters L* , that is, for any contour 7 e F 
in the domain, 7 0(7), where 0(7) is a number of absolute value one, and for 
any two contours 71 , 72 

a (71 72) = a(7i)a(7 2 ). 
The Hardy space H°°(X + , a) consists of holomorphic multivalued functions / uni- 
formly bounded in the domain C\E such that 

f(jz) = a(j)f(z), 

and H 2 (X + , a) is the closure of H°°(X + , a) with respect to the norm (5.9). 

Note that the absolute value of a function from H 2 (X + ,a) is single valued and 
a fixes, actually the ramification of the argument of the function. 

Example. As it was mentioned, the function b is not single valued but \b(z)\ is a 
single valued function. We define the character /j, £ T* by 

Let 7j be the contour, that starts at infinity (or any other real point bigger than 
1), go in the upper half-plane to the gap (aj,bj) and then go back in the lower 
half-plane to the initial point. Assuming that bo < ... < aj < bj < < ... < &o, 

we have fJ,(~/j) = e~ 27Tl 9+ 1 , equivalently u([bj,ao]) = g g+7 J ■ 

Remark. Note that the system of the above contours jj is a generator of the free 
group T*(E). In other words a character a is uniquely defined by the vector 

[a(7i), a( 72 ), a( lg )]£V. 

This sets an one-to-one correspondence between T*(E) and T 9 . 

Proposition 5.4. Using the above definitions we get the orthogonal decomposition 

H 2 (X + ) = {1} H 2 (X + ) = {1} 8 bH 2 {X +ll i- 1 ). (5.12) 

Now we can iterate (5.12). Let k a be the reproducing kernel of H 2 (X + , a) with 
respect to infinity, that is, the vector from H 2 (X +7 a), which is uniquely defined 
by the condition 

(f,k a ) =/(oo), Vf£H 2 (X + ,a). 

Then 

h 2 (x+) = {1} e {bk^ 1 } e b 2 H 2 (x+, M - 2 ), 

and so on... 

Theorem 5.5. Let a £T*. The system 

K}„ ez+ , e«:=b n K<*»-\ ^-" : =^_^ (5.13) 

forms an orthonormal basis of H 2 (X + , a) . The same system with n £ 7L is an 
orthonormal basis in L 2 (dX + ). Moreover, the multiplication operator by z is a 
periodic Jacobi matrix with spectrum E. 

Theorem 5.5 indicates a special importance of the functions k a . They are very 
well studied [6]. First of all, they have analytic continuation (as multivalued func- 
tions) on the whole X, so we can write k a (Z). 
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Theorem 5.6. For every a € T* the reproducing kernel k a {Z) has on X exactly g 
simple poles that do not depend on a and g simple zeros. The divisor D + = J2j=i Zj 
of zeros 

k a (Zj)=0 (5.14) 

with the divisor of poles form the divisor 

div(k a ) =D+-D C (5.15) 

that belongs to D(E), moreover (5.15) sets an one-to-one correspondence between 
D{E) andT*{E). 

The functions k a possess different representations, in particular, in terms of 
theta-Tunctions [9], and the map D i— > a can be written explicitly in terms of 
abelian integrals (the Abel map). 

Summary. The three objects J(E), D(E) and T*{E) are equivalent. Both maps 
T*{E) — > D(E) and T*(E) — ► J(E) can be defined in terms of the reproducing 
kernels of the spaces H 2 {X +1 a) 1 a £ T*(E). The first one is given by (5.15). It 
associates to the given k a (Z) the sets of its zeros and poles (the poles are fixed and 
the zeros vary with a). The matrix J (a) e J(E) is defined as the matrix of the 
multiplication operator by z(Z) with respect to the basis (5.13): 

z(Z)e»(Z)=p«e«_ 1 (Z) + q«e»(Z)+p« +1 e» +1 (Z), Z e X, s£ Z. (5.16) 

It's really easy to see that J(a) is periodic: just recall that b 9+1 is single valued, 
that is, [i 9+1 — 1, and therefore the spaces H 2 {X + , a) and H 2 (X + , a^^ 9+1 ^) (and 
their reproducing kernels) coincide. 

Now we can go back to the Renormalization Equation. Note that 7r acts naturally 
on T{E): 

tt 7 = {tt(Z), Ze~f}e T(E), for 7 e T(E). 
The map tt* : T*(E) -» T*(E) is defined by duality: 

(ir*a)(j) = a(n-f). (5.17) 

Theorem 5.7. Let T , T _1 : [—1, 1] — > [—1, 1], be an expanding polynomial. Let J 
be a periodic Jacobi matrix with spectrum E C [—1,1], and therefore there exists 
a polynomial U such that E = C/ _1 [— 1,1] and a character a <G T*(E) such that 
J = J (a). Then 

J := J(n*a) (5.18) 

is the periodic Jacobi matrix with spectrum E = C/ _1 [— 1,1], U := U oT, that 
satisfies the Renormalization Equation (2.1). 

Proof. First we note, that for the operator multiplication by z(Z) in L 2 (dX + ), the 
operator multiplication by z(Z) in L 2 (dX + ) 1 the spectral parameter z and the 
isometry 

(vf)(Z) = f(n(Z)), v : L 2 (dX + ) - L 2 {dX + ), 

we have 

/ ; \ (7) \W)(Z)\ 2 d»= [. \\ E 7 \ (7] )\f(Z)\ 2 du;. (5.19) 
Jdx + zo-z(Z) Jdx +\ v (?)Lz ) 
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It is evident, that 

_ T'(z )/d 



1 _J_ 

d Zn - 



T(y)=x z o-y T(z )-x- 

Thus 

v*(z - z{Z)Y x v = (T'(z )/d)(T(z ) - ~z{Z))-\ (5.20) 
It remains to show that v transforms the basis vector 

-a = ln K &ii- n 

into 

<d a = b nd K^' a ^~ nd = (b n o n)K 7r *( &ii ~ n) . 

Or, what is the same, that K 5 o tt = K** a for all a e r*(S). Note that both 
functions are of norm one in the same space H 2 (X + ,n*a), in particular, they have 
the same character of automorphity n*a £ T*(E). Note, finally, that the divisor 

div(fc" O TT) = TT- 1 {D+) - TT^iDc), 

where div(fc a ) = D + — Dc, belongs to D(E), therefore k a o n is the reproducing 
kernel and the theorem is proved. □ 

To find all other solutions of (2.1) let us look a bit more carefully at the above 
proof. 

Note that the same identity (5.19) holds for any isomctry v of the form 

vf = v e f = 9(f ott), 

where 9 is a unimodular (\6\ = 1) function on dX + . 

Concerning the second part of the proof, let us mention that the set of critical 
points of U splits in two sets: 

{c u }=T- 1 {c }u{c T }. 

Correspondingly, 

E( C ^ = E E {c u ) kJ +'E{c T ) j , 

k n { Cu) k ,i={Cuh 

and the divisor of k a o n consists of two parts, that one that depends on a 
and that part that corresponds to the critical points of the polynomial T 
since 

d-l d-l 

D = div(fc« o tt) - n-^D) + E(^t), - n-^Dc) - E^" 

Thus we can fix an arbitrary system of points {Z c j}^~l such that z(Z c j) belongs 
to the same gap in the spectrum E as the critical point {cr)j- If 9 is the canonical 
product on X with the divisor 

d-l d-l 

div((9) = E z c,i-E( c 'r)i, 

3=1 3=1 
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then 0k a o tt is the reproducing kernel simultaneously for all a G T*(E). But to 
make 9 unimodular (zeros and poles are symmetric) our choice is restricted just to 
Z c ,j — {Ct)o or Z c ,j — t(Ct)j- Note that T(Cr)j — (Ct)j is the devisor of the 
Complex Green function &(c r )j - In this way we arrive at 

Theorem 5.8. For an expanding polynomial T, and a periodic Jacobi matrix 
J = J (a), a E T*(E) as in Theorem 5.7 there exist 2 d ~ 1 solutions of the Renor- 
malization Equation (2.1). Denote by fJ-( CT )- the character generated by the Green's 
function b( CT }., &( Cr )j °7 = l 1 (c T ) j (l)b(c T ) r Then these solutions are of the form 

J :=J(r, S ir*a) , ns := U^^, (5-21) 

as before 

6 = {^(ct),}, 5 {ct)] = ±1. 
Proof. We define the isometry 

(vf)(Z) = (Ub^'^\z)^ f(n(Z)) 

and then repeat the arguments of the proof of Theorem 5.7. □ 

Concluding this section note that the central part in the proof of Theorem 1.2 
(the claim that the limit matrix has a form of the orthogonal sum) also can be 
reduced to an another well known fact from the theory of Hardy spaces on Riemann 
surfaces. Namely to the statement that H 2 is trivial, i.e., 

H 2 (C\E) = {const} (5.22) 

in a domain of the form C\E, where the Lebesgue measure of \E\ = 0. 

An alternative proof of (4.1). Since we can start with an arbitrary J, we start with 
a periodic matrix related to a certain H 2 (X + ), e.g., with the matrix with constant 
coefficients J = E = [—1, 1]. Then, under inverse iterations of the polynomial 

T according to Theorem 5.7, we will get spaces of the same nature (i.e., the character 
is trivial, equals one on every contour). Let J n be the matrix with spectrum E n = 
(T™) _1 [— 1, 1]. We have 

z(Z)e 1 (Z,n) = p(n)ie (Z,n) +q(n)iei(Z) +p 2 e 2 (Z,n), Z e X n , 

here n is related to the number of iterations and the position of the element of the 
matrix is fixed. Recall that eo(Z,n) = 1 (the initial basic vector, see (5.12)) and 
we have, putting Z — oo, 

p(n)i = (zb n )(co)K^™ (oo,n). 

(z6„)(oo) is the so called capacity of E n , if it goes to zero even better, in fact it 
does not, but in any case it is uniformly bounded. Then, assuming that if M ™ (oo, n) 
does not go to zero, by compactness arguments, we can find a subsequence 

(b nj )(z)K^ (z, nj ) 

that converges pointwise in the domain to a non-trivial holomorphic function from 
H 2 (C \ E), E — limi^n = Julia(T), that equals zero at oo. But this contradicts to 
(5.22). Thusp(n)i -> 0. □ 
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Note that this proof is valid for expanding polynomials (we do not require that T 
is sufficiently hyperbolic). Note also the flip in notations of the matrices' elements: 
in this section a basis of holomorphic functions substitutes the standard polynomial 
basis (instead of the multiplicity of the pole at infinity we enlarge the multiplicity 
of zero). That is, p\ in this section is the same as po in Section 4 (we are just unable 
to enumerate the elements, related to holomorphic functions, by negative integers). 

6. Concluding remarks 

Our concluding remarks concern basically other solutions of the Renormalization 
Equation. 

6.1. The duality S i— > — S. In Theorem 1.1 we proved contractibility of only one 
of the solutions of the renormalization equation corresponding to S = 5- , but it 
means that at least one more solution has the same property. 

Theorem 6.1. The dual solution of the Renormalization Equation J(J,—5), pos- 
sesses the contractibility property simultaneously with J( J,6). 

It deals with the following universal involution acting on Jacobi matrices 

J — > J T := U T JU T , where U T \l) = |1 - I). (6.1) 

Obviously VU T = U T S 1 ~ d V. Thus, having J as a solution of the renormalization 
equation corresponding to J we have simultaneously that S d ~ 1 J T S 1 ~ d solves the 
equation with the initial J T . The following lemma describes which branch corre- 
sponds to which in this case. 

Lemma 6.2. Let J = J( J, S) then 

S d - 1 J T S 1 - d = J(J T ,S). (6.2) 

Proof. We give a proof using the language of Sect. 5, so formally we prove the 
claim only for periodic matrices. 

Note that the involution (6.1) is strongly related to the standard involution r 
(5.2) on X. Indeed, the function K(tZ, a) has the divisor 

tD + - tD c = (tD+ - D c ) - (tD c - D c ), 

that is, 

K(Z,(3) 

K{TZ > a) = b cl{ z)... bcg{ zy 

and (3 = va^ 1 , where v = \i Cl . . . fi Cg . Due to this remark and the property z(tZ) = 
z(Z) we have 

(J(a)) T = J^iia- 1 ). (6.3) 
Now we apply (6.3) to prove (6.2). Let J T = J (a) with a G T*(X + ). Or, in 
other words, J = J(pi)a^ 1 ). Then by (5.21) 

d_1 i (1+5 ) 

J{J,5) = Jins^ifiva- 1 )), r] S := Y[^ CT)j ^ • 

But 77* ft = p d and 7r*(/>) = vq^ (just to look at the characters of the corresponding 
Blaschke products). Thus, having in mind that ris^s — VS+, we obtain 

j(J,s) = j(Ac*V(0)- 
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Using again (6.3) we get 

(J(J,S)) T = JO^-V^)) = 5 1 - d J(r ? _ 57 r*(a))5 d - 1 , 
and the lemma and Theorem 6.1 are proved. □ 

Having two different contractive branches of solutions of the renormalization 
equation, following [8] , to an arbitrary sequence 

e = {eo,£i,£2 • • • }, £j ; = S±. 
we can associate a limit periodic matrix J with the spectrum on Julia(T). For a 
fixed sufficiently hyperbolic polynomial T, we define J as the limit of 

J n := J(7] eo TT*7] ei ...TT'lfc^J. (6.4) 

6.2. Other solutions of the Renormalization Equation and the Ruelle op- 
erators. We conjecture that actually all branches of solutions of the renormaliza- 
tion equation are contractions for sufficiently hyperbolic T. At least the previous 
remark looks as a quite strong indication in this direction: considering, instead 
of initial T, T 2 = T o T or its bigger powers, we get, as in (6.4), several S's, 
Vs = Veo^Ve! ■ ■ ■ 7r *?7e„-i7 possessing the contractibility property with respect to 
the polynomial T n and different from 5± (related to {■K*) n ). 

Similarly to (4.4), (4.5) we formulate 

Conjecture 6.3. Let T(z) be an expanding polynomial and let T'(z) = A\(z)A2(z) 
be an arbitrary (polynomial) factorization of the derivative. Denote by <J\,2, the 
(nonnegative) eigen-measures, corresponding to the Ruelle operators 

T(y)=x iyy > 

i.e., L* A .<Ji — piOi. Finally let Ji,2 be the one-sided Jacobi matrices associated with 
(Ti,2- Then the block matrix J = J_ ® J + with J_ = J\ and J + = J 2 is limit 
periodic. 

Note that by the same reason as above the conjecture holds true, say, for T 2 (z) 
and - T'(z), A 2 = T'(T(z)). 

6.3. Shift transformations with the Lipschitz property. We say that the 
direction r\ £ T* has the Lipschitz property with a constant C(r}) if for all a, (i e T* 

\\J( V a)-J( V P)\\<C(v)\\J(a)-J((3)l (6-6) 
Then, one can get the contractibility of the map rjir* in two steps: 

\\J( V n*a) J( V n*P)\\ < C( V )\\ J(n*a) - J{„*0)\\ 

<C{ V )x\\J{a)-J0)\\. 

Note, that in fact the situation is a bit more involved because we should be able 
to compare Jacobi matrices with different spectral sets, for example, when Ei — 
T~ 1 E i , Ei 7^ Ei. But we just wanted to indicate the general idea, in particular, for 
directions rjs of the form (5.21) such a comparison is possible. Of course, for our 
goal the constant C{rf) should be uniformly bounded when we increase the level of 
sufficient hypcrbolicity of T making x smaller. 
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However the key point of this remark (this way of proof) is that, actually, we 
do not need to constrain ourselves by the form of the vector rj. Combining a "Lips- 
chitz" shift by r\ (the direction is restricted just by this property) with a sufficiently 
contractive pull-back ir* we arrive at an iterative process that produces a limit peri- 
odic Jacobi matrix with the spectrum on the same Julia(T). In the next subsection 
we give examples of directions with the required property, see Corollary 6.6. 

We do not have a proof of the Lipschitz property of rjs 's, but there is a good 
chance to generalize the result of the next subsection in a way that at least some 
of the directions rjs will be also available. 

Finally, we would be very interested to know, whether there is in general a 
relation between the form of the "weight" vector 77 and the corresponding weights 
of the Ruelle operators (if any exists). 



6.4. Quadratic polynomials and the Lipschitz property of the Darboux 
transform. Consider the simplest special case T(z) = p(z 2 — 1) + 1, p > 2. Note 
that the spectral set E = T~ X E is symmetric, moreover the matrix related to 
H 2 {-K*a) has zero main diagonal (as well as a one-sided matrix related to a sym- 
metric measure). Now we introduce a decomposition of H 2 (ir*a) which is very 
similar to the standard decomposition into even and odd functions. 

We define the two-dimensional vector-function representation of / € H 2 (n*a) 



1 

71 



'f(Zi(z)y 


1 — ► 


'gi(z) 


J(Z2(Z))_ 




MZ). 



(6.8) 



where 



'gi(z) 


1 


MZ). 


_ 2 



/(Z 1 (Z)) + /(Z 2 (Z)) 

f{z 1 Cz))-f{z 2 {z)) 

the first component, in a sense, is even and the second is odd. To be more precise, 
let us describe analytical properties of this object in dX + . 
Note that due to 



/ l/l 2 ^- I. \ 

JdX + JdX + L 



E 

tt(Z)=Z 



\f\ 2 dQ 



metrically it is of L 2 with respect to w, moreover the transformation is norm- 
preserved. 

It is evident that the function 171 belongs to H 2 (X +7 a). Consider the second 
function. Note that the critical points of T are zero and infinity. For a small circle 
7 around the point T(0) = — p + 1 we have 32 7 = —92 and the same property for 
a contour 7 that surrounds infinity. Let us introduce 

A 2 := b T(0) b. 

Note that for the above contours we have A o 7 = — A. We are going to represent 
g 2 in the form g 2 = Ag 2 and to claim that g 2 has nice automorphic properties in 
X + . Let us note that 

T(0) l2 z-T(0) 



b Z -^ 



°T(0) 



A 2 
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is an outer function in the domain C\E ~ X+. So, the square root of this function 
is well defined. We put 



and denote by fj the character generated by 0, (^07 = 77(7)0. Thus (6.9) reduces 
the ramification of the function A to the function (p, which is well defined in the 
domain, and to the elementary function y/z — T(0). 

Theorem 6.4. The transformation / ^ ffi © 32 given by (6.8) is o unitary map 
from H 2 (n*a) to H 2 (a) H 2 (afj). Moreover with respect to this representation 



zf 

and 



"0 cf 




51" 


4> 




92_ 



(6.10) 



v+f^f®0, feH 2 (a), 
where the isometry v+ : H (a) — > H 2 (ir*a) is defined by (5.11). 
Proof. By the definition of A we have 

,g 2 = Ag 2 , where g 2 € H 2 (aff). (6.11) 

Further, since 



21,2 = ±1 

we have, say for the second component, 



T(0) 



= ^.9i- (6.12) 



1 (z/)(Z(ZQ)-(z/)(Z(Z 2 )) _ z-T(0)f(Z(Z 1 )) + f(Z(Z 2 )) 
A 2 y pA 2 2 

Since on the boundary of the domain 

(the second expression is positive on dX + ) we have 

(f>A 2 = ~§ on E. (6.13) 

Using this relation, similarly to (6.12), we prove the identity of the first components 
in (6.10). 

□ 

Theorem 6.5. The multiplication operator <j> : L 2 (dX + ) — > L 2 (dX + ) with respect 
to the basis systems (5.13) related to a and f\a, respectively, is a two diagonal 
matrix $. Moreover, 

M = M , M > = J ^- r (°) , (6.14) 
P P 

In o£/ier words, the transformation J (a) 1— > J{rja) is the Darboux transform. 
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Proof. First of all <f> is a character-automorphic function with the character fj with 
a unique pole at infinity (b(f) is an outer function). Therefore the multiplication 
operator acts from bH 2 (ap~ 1 ) to H 2 (rja). Therefore, the operator $ has only one 
non-trivial diagonal above the main diagonal. The adjoint operator has the symbol 
4>. According to (6.13) it has holomorphic continuation from the boundary inside 
the domain. Thus is a lower triangular matrix. Combining these two facts we 
get that <3? has only two non-trivial diagonals. Then, just comparing symbols of 
operators on the left and right parts of (6.14), we prove these identities. □ 

Corollary 6.6. Let Ji, 2 be periodic Jacobi matrices with the spectrum on [—1,1]. 
Let Darb(Ji j 2,p) be their Darboux transforms. Then 

||Daxb(ii,p)-Darb(J 2 ,p)|| < C(p) || J x - J 2 1| . (6.15) 

Proof. For the given Ji ; 2 we define Ji j2 via the quadratic polynomial T(z) — p(z 2 — 
1) + 1. Being decomposed into even and odd indexed subspaces they are of the form 





$1,2 



*1,2 





(6.16) 

2 1| , we have 
(6.17) 

with n(p) = -32, C is an absolute constant. Using (6.14) wc get (6.17) with 

c( P ) = " ' a 



Jl,2 = 

Due to the main theorem, that gives the uniform estimate for || J\ 

||*i -* 2 || <«(p)||Ji- M, 

if ic on o Kon 

p-2 

P=2 

7. Appendix 

Here we recall some basic facts on two-sided Jacobi matrices. Let J define a 
bounded selfadjoint operator on Z 2 (Z). The resolvent matrix-function is defined by 
the relation 

-(0\(J-z)-i\0) (0\(J - z)-i\iy 

_(l\(J - z)-i\0) (l\(J-z)-i\l)_ 
This matrix-function has an integral representation 

da(x) 



W{z) = W(z,J) = 



(7.1) 



W(z) 



(7.2) 



with 2x2 matrix-measure having a compact support on R. J is unitary equivalent 
to the multiplication operator by an independent variable on 

'fo(x) 



K = \f = 



{'-ni™<~>} 



(7.3) 



moreover, under this unitary mapping from I 2 — > L 2 we have 

ID ^ 



(7.4) 



Let r-{z) = r-(z, 0), r+(z) = r + (z, 1) be resolvent functions of J_ = J_(0) and 
J + = J+(l) respectively (see (2.2)), and a± be the corresponding (scalar!) spectral 
measures. Then 

rZ 1 (z) pi 
Pi r+ 1 (z)_ 



W(z) 



(7.5) 
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Recall that according to our notation piQd is the orthonormal polynomial of the 
degree d — 1 for the measure a+ and piRd is the related polynomial of the second 
kind: 

x PiQd{x) -piQd(z) 



PiRd(z) = / da+(x) 



In this notations 
moreover 

Lemma 7.1. Let 

Then 



\d) i ^ £ d {x) :-- 



-p\R d {x) 
PiQd{x) 





PiRd(x) 



1 -p?i2d(x)' 
piQd(x) 



da{x) 



£ d (x) - £ d (z) 





p\Rd{x) 



(7.6) 
(7.7) 
(7.8) 

(7.9) 



(0\(J - z)-'\0) (0\(J-z)-i\d) 
(d\(J - z)-i\0) (d\(J - z)~i\d) 



T*(z)W{z)T(z) + H(z)G(z). (7.10) 



Proof. This is a standard trick from the theory of orthogonal polynomials. Due to 
the unitary mapping onto L^, equivalently we have to calculate 

T{x)*da{x)T{x) 



I 



x — z 



(7.11) 



Therefore, using orthogonality and (7.8), we continue 



F{x)* - Hz)* 



=F{z)* j da(x) 



da{x)Hx) + T(z 



T(x) - T{z) 



x — z 

=T*{z)g(z)+T*(z)W{z)T(z). 



T(z) 



7 
7 



da(x)J r (x) 



da{x) 



Hz) 



(7.12) 



□ 



Corollary 7.2. Combining (7.5) with (7.10) we get (2.23) from 

\(Q\(J-T(z))^\Q) (0\(J-T(z))-i\l) 
(l\(J-T(z))-i\0) (l\(J-T(z))-i\l\ 



(okj-^io) (o\(j- z yi\d) 

{d\(J - z)-i|0) (d\(J - z yi\d) 
which is a part of the Renormalization Equation. 
Proof. A straightforward computation. 



T'(z)/d 



□ 
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